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Subthermal switching with nanomechanical relays 

Dincer Unluer and Avik W. Ghosh 


Abstract —We present a physical model for electronic switching 
in cantilever based nano-electro-mechanical field effect transis¬ 
tors, focusing on the steepness of its switching curve. We find 
that the subthreshold swing of the voltage transfer characteristic 
is governed by two separate considerations - the ability of the 
charges to correlate together through dipolar interactions and 
amplify the active torque, versus the active pull-in forces that 
drive an abrupt phase transition and close the air gap between the 
tip of the cantilever and the drain. For small sized relays, dipolar 
and short-range Van Der Waals ‘sticking’ forces dominate, while 
for longer cantilevers the capacitive energy acquires a major 
role. The individual pull-in and pull-out phases demonstrate a 
remarkably low subthreshold swing driven by the capacitive 
forces, sharpened further by dipolar correlation. The sharp 
switching, however, comes at the expense of strong hysteresis as 
the metastable and stable states interchange along the forward 
and reverse phases of the voltage scan. 

Index Terms —Nano-electro-mechanical systems (NEMS), 
nano-electro-mechanical field effect transistor (NEMFET), sub¬ 
threshold swing, van der Walls forces, hysteresis. 

1. Introduction 

D igital electronics is fundamentally constrained by the 
cost of binary switching QAV. The number of charges 
N = Q/q is quite large (between 10-100 in the past 10 years) 
for a n-MOS transistor, and ramps up rapidly to a few thousand 
for a fan-out 4 device with interconnect, because of the need 
to charge up the entire interconnect capacitance and drive an 
analogous transistor downstream (D HI AV on the other 
hand is constrained by the Boltzmann limit, and corresponds to 
a switching slope of ksT In lO/g' ^ 60 mV/decade. Two paths 
to bypassing these constraints while maintaining an acceptable 
reliability are (i) correlated switches as in nano-magnetic logic 
ca, where several units, in this case spins, can offer solidarity 
against thermal fluctuations, and yet limit the energy cost 
by correlating together through their internal coupling fields; 
and (ii) phase-transition switches that operate either far from 
equilibrium or near a phase transition point, so that a small 
voltage precipitates a large change, such as through an internal 
voltage amplification. 

In order to lower the power supply voltage Vdd and thereby 
cut static and dynamic power dissipation, we ultimately need 
devices with lower sub-threshold swings and steeper gate volt¬ 
age transfer curves beyond what the Boltzmann limit posits. 
This has led to the exploration of various unconventional 
switching devices such as tunnel field effect transistors (H, im¬ 
pact ionization devices ISl, ferroelectric transition switches (6) 
and chiral tunneling devices based on graphene pn-junctions 
Q, all of which purport to achieve swings lower than the 
Boltzmann limit. In this regard, mechanical switches or relays 
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are noteworthy, as they can swing between near zero off- 
currents and high on-currents at very low Vdd values IH. Such 
nano-electro-mechanical field effect transistors (NEMFETs) 
are relevant as low power switches El, especially when 
integrated into hybrid CMOS circuits to minimize leakage 
current from the power supply supply in the OFF state CQl. 
Indeed, ultra-low subthreshold swings have been demonstrated 
using suspended-gates ifTTIl . and micro-relays 13, (121, along 
with complimentary relays for logic applications with larger 
voltage margins and reliability ca. There now are multiple 
designs including the printable relay ca and the curved relay 
by IBM Ha. Relays are also known to operate as active 
elements in turning on voltage gated ion channels in axonal 
systems. In addition, it is noteworthy that the voltage gating of 
sodium channels involving a mechanical relay was measured 
many decades ago at a subthreshold swing of 20 mV/dec ca, 

ITTl . ITsl . 

The action of a NEMFET is significantly different from 
a conventional electrostatically driven field effect transistor, 
and needs new physics-based predictive models. A number of 
empirical models have already been proposed based on ID 
electrostatics lEl, OSl, nonlinear oscillator equations (201, 
fm . (221, pull-in instability (23 . (241, deformation char¬ 
acteristics of microstructures under electrostatic loads (23, 
and techniques from chaos and bifurcation theory in finite- 
elements (261, (23, (2S1- Our purpose in this paper is to 
explain the electron transport physics in a relay switch, starting 
with the microscopic Landauer equation that also elegantly 
explains the limits on the sub threshold swing. We argue that 
a dipole coupled, cantilever based relay can function as both 
a correlated switch and a phase transition switch. 

II. Device Structure and transport equations 

Nanoelectromechnanical switches exist in various fiavors - 
movable gates or movable channels. Movable gate switches are 
driven by the capacitive energy CVqI2, where the capacitance 

C = a(^ox/^ox + IS divided serially between 

a fixed oxide and a variable air gap. Upon application of a 
gate bias, a force F = —dUjdVG alters the air gap from its 
equilibrium value to, counterbalanced by a spring which exerts 
a restoring force —k{tair—to) (23 • At a critical distance where 
the electrostatic pull-in overcomes the restoring force, there is 
an abrupt closing of the air gap, precipitating a sudden rise in 
capacitance. Since the pull-in acts on a capacitive structure, 
one can play with material properties of the gate oxide and 
shape its potential landscape in a creative way in order to 
liberate the movable gate easily (30l . Since the channel itself 
is fixed, its OFF current is limited by electrostatic depletion, 
and for a small channel suffers from the usual source-to-drain 
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Fig. 1. (Left) Schematic picture of cantilever showing all the forces operating 
on it. (Right) Associated energy terms, spelled out in Eq. 


tunneling, although the removal of the movable gate reduces 
the gate oxide leakage significantly. 

In a movable channel transistor, the idea is to physically 
separate the channel away from the drain so that in addition 
to electrostatic depletion, the OFF current can be reduced sub¬ 
stantially with an added quantum mechanical tunneling barrier. 
In fact, it is straightforward to show that their electrostatic 
and conformational modes and thus their subthreshold swings 
act in parallel, S~^ ^ tSTl . assuming a clear 

separation of electronic and conformational time constants. 
We have shown in the past ED that the conformational mode, 
besides aiding the electrostatic turn-off, allows an intrinsically 
high transconductance, and thus a low subthreshold swing 
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where (7 is a correlation function that we will describe later. 
The cantilever FET can function as both a correlated switch 
and a phase transition switch. The presence of the dipole /i 
means that the subthreshold swing can be reduced by a factor 
Qtox! octox/L for longer cantilevers, exploiting the property 
that many charges are physically moving with the cantilever 
for the price of one. This is, in fact, the suggested mechanism 
behind the observed sharp switching in voltage-gated sodium 
channels ca. A separate mechanism not analyzed in ED 
arises from an additional ‘pull-in’ that sucks in the end of 
the cantilever once it is physically proximal with the drain. 
The short-ranged Van der Waals force and the longer ranged 
capacitive forces cause an abrupt shrinkage of the air gap and 
an exponential rise in tunneling current, reducing the sub¬ 
threshold swing further by a factor proportional to dtairldO, 
which vanishes near the contact point, giving a near vanishing 
subthreshold swing. 

The Van der Waals pull-in helps subvert thermal fiuctua- 
tions by stabilizing the attractive forces. Resetting the device 
however, requires circumventing adhesion, which creates a 
sizable hysteresis loop. The switch thus has extremely sharp 
transitions along the walls of the hysteresis curve, but its loop- 
averaged subthreshold swing is quite large and problematic. 
Fixing this requires creative material and geometry design. 
Our purpose in this paper is to develop the mathematical 
model and the physical picture behind the entire cantilever 
dynamics in its ON, OFF and transitory states, the geometry 
of the hysteresis and its dependence on material parameters as 
well as the voltage scan rate. 

The electron current is determined by the degree of quantum 
mechanical tunneling across the gap between the cantilever 


and the drain. We can describe that using the Landauer 
equation ES. written in an equivalent but less familiar ED 
form as a convolution between the zero-temperature current 
and a thermal broadening function Ft{E) = —df{E)/dE 

dEM{E)(f{E, Vg)^ ® Et{E, Vg) (2) 

where M{E) is the mode count in the cantilever, / is the 
equilibrium Fermi-Dirac distribution, with the Vq dependence 
in Ft arising from the location of the Fermi energy Ef 
which shifts with Vg scaled by a gate transfer factor ac- 
Pi ^2 are the bias separated electrochemical potentials in the 
source and drain contacts, and T{E,Vg) is the mode and 
conformation averaged transmission function controlled by 
the gate voltage-dependent orientation. (...) represents the 
thermodynamic average of the cantilever orientation 0 over 
various angles, weighted by its probability distribution P 

, , f Mt(E,e)P{t),Va) 

If we assume the mechanical modes in the cantilever are at 
equilibrium, (i.e., the potential is ramped in time much slower 
than the relaxation scales in the lattice), P is simply the 
Boltzmann solution oc exp [—U{0,VG)/kBT], even though the 
electrons flowing are clearly in non-equilibrium. In general, 
the cantilever’s angular distribution satisfies a thermodynamic 
equation that needs to be derived properly from a stochastic 
Langevin equation in the presence of a fast ramp speed of 
the potential. A fiexible cantilever clamped on one end and 
free on the other satisfies a fourth degree partial differential 
equation of motion. For an over-damped, rigid cantilever with 
a single defining angle 0, the equation can be simplified down 
to second degree and consequently P can be shown to satisfy 
the Fokker-Planck equation 


dP d ( dU{d,VG) \ 
dt 89 y 89 j 




d‘^p 


(4) 


where 7 is the angular damping constant with units J-s. We 
assume here that all phonons emitted in the cantilever are 
promptly reabsorbed by a conductive thermal coating coupled 
to the source, so we can ignore any self-heating and phonon 
bottleneck effects in the cantilever itself. 

The above equations carry all the necessary ingredients 
for our NEMFET operation. The normal subthreshold swing 
of kBTlnl^/acq arises from the gate dependent shift in 
Fermi energy relative to the mode placement M for a fixed 
orientation, when f/ is a delta function around a given 0, and 
P is simply equilibrium Boltzmann. Further reduction in sub¬ 
threshold swing will arise from the added gate dependence of 
the conformational modes, hidden in U{0 ,Vg)- The hysteresis 
will arise from a transition out of metastable states with a non¬ 
equilibrium distribution P in presence of a time-dependent 
sweep in gate voltage, Veit). Finally, a negative capacitance 
for a ferroelectric oxide EOl will be captured by the gate 
voltage dependence of the oxide dielectric constant, eoxiVc) 
sitting in U, which can introduce its own dynamics near the 
ferroelectric phase transition point. 
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III. Operation of the relay 


The physics of the relay is captured by its transmission 
and its conformational potential profile. The transmission 
is limited by electron tunneling into the drain from the 
edge of the cantilever through the vacuum barrier of width 
tair{^) = h — LsinO, h being the bare height above the 
drain before bending and L being the length of the cantilever 
(Fig. 0- The transmission, using a WKB approximation, is 
T{E,0) (X exp[—2Ktair{0)]. The main physics of the NEM- 
FET dynamics is captured by its potential U{0,Vg), which 
can be separated into a flexural potential Uq representing 
the tendency of the cantilever to snap back to its horizontal 
orientation, a gate driven field £ which couples with the dipole 
moment jl along with capacitive coupling with the back gate 
to bend the cantilever, and finally a Van Der Waals ‘pull-in’ 
potential Uvdw that causes the cantilever to ‘stick’ to the 
drain when close enough (Fig. 

U[0,Vg] = 

Uo{0) ^ 

UvDw{^) = 

UdiplO^VG] = 

UcaplO^VG] = 


Uo{0) + UvDw{0) 
Udip[0,VG]EUcap[0^VG] 

lu. 


2 

4e 


Bendn‘2 
0 ^ 
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tair (^) 
-il-£ = -m| 
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tair (^) / 
Vg sin 0 


^oxj^ox i^air 
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We assume a simple parabolic potential for the bending 
dynamics to focus on the basic physics. Higher nonlinear 
terms would capture the full import of the fiexural modes 
El, although they do not provide us with immediate insights 
ED. It is worth emphasizing that the dipolar torque has a 
distinct advantage over capacitive forces, as it exploits the 
vector directionality of the gate field. Even if the oxide is 
thick enough to normally produce short-channel effects, the 
gate field is orthogonal to the drain field and thus creates a 
larger torque than the drain field (especially in the presence of 
dual gates), which makes its effect more pronounced. As the 
device sizes grow towards the micron regime, the increasing 
width (W) starts to overcome sticking, so the effect of the 
capacitive forces starts to overcome the effect of the Van Der 
Waals and dipolar forces. 


TABLE I 

Simulation Parameters eor Model 


L 

5 nm 

keT 

0.025 eV 

h 

2.5 nm 

jjBend 

2 

W 

0.5 nm 

P 

0.15 eA 

i-air 

1.5 nm 

^Relay 

3.98 eV 

^air 

1 

^MetalPlate 

0.67 eV 

tox 

1 nm 

^ Relay 

3.47 A 

^ox 

3.9 

^ MetalPlate 

3.144 A 



Fig. 2. (A) The parabolic potential of the cantilever with one local minimum 
resulting from the combination of the stiffness and the external field terms. 
The local minimum shifts with the gate bias. (B) The VDW potential resulting 
from the gate electrode with a stationary local minimum. (C) A microscopic 
view of the parabolic potential varying with gate bias, in addition to a localized 
VDW potential from the contact or the capacitive destablization, creates a net 
profile that has one gate-dependent local minimum and one relatively fixed 
global minimum. 


Fig. shows the fiexural parabolic potential, progressively 
shifted by the external electric field, causing a change in 
average angular orientation. Fig. 1^ shows the separate Van 
der Waals potential that acts as a pull-in when the relay 
comes close enough to the drain. In combination of these two 
effects (Fig. |^), there is a sweet spot where the potential 
profile between the metastable oscillator minimum and the 
stable Van der Waals minimum reaches a point of infiection, 
whereupon the cantilever angular coordinate makes an abrupt 
phase transition and locks onto the drain (in actuality this 
happens a bit earlier, once the hill separating the Van der Waals 
and parabolic minima reduces to around ^ The pull-in 

has an important role to play. In its absence, the cantilever 
would be either too rigid to displace with an external field, 
or too floppy to resist thermal averages. In the absence of 
any other barrier, the thermal average over the conformational 
potential would tend to restore the system to its initial OFF 
state. This is a well known problem with a single well, analog 
computing - i.e., while it is energetically very pliable, it is 
highly corruptible with noise. In our geometry, we have an 
initial low energy analog variation combined with dipolar 
coupling, followed by a digital component at a later stage 
of the game, once the system abruptly falls into its ON state 
configuration in the well. In other words, it provides only a 
weak barrier to forward motion from the metastable to a stable 
state but a strong barrier for a return to origin. 

The barrier asymmetry is at the heart of the NEMFET 
operation, but relies on the intrinsic non-equilibrium nature of 
the dynamics. Left to itself long enough, the cantilever would 
always find a way to spontaneously jump out thermally from 
its metastable state to the stable state, i.e., it would bend over 
and stick to the drain given adequate time even in absence of 
a gate field assisting it. To define the transition between states, 
we will need to operate the voltage scanning faster than this 
spontaneous transition rate, so that the cantilever coordinate 
stays pinned to the metastable state until the voltage eliminates 
the barrier to the stable ground state. Our transport model 
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Fig. 3. Pull-in (shown by blue curves) and pull-out (shown by the red curves) 
mechanism rising from the NEMFET potential at positive and negative gate- 
biases. (Right) The pull-in occurs when the positive gate-bias lowers the local 
maximum (barrier) between the two local minimums allowing the cantilever 
to move right. (Left) The negative voltage shifts the local minimum from the 
VDW forces to higher energies, thus allowing the cantilever to shift left to 
the local minimum from the parabola. 


must also probe this intermediate regime, and a simple steady 
state model will not capture the hysteresis associated with the 
metastable state. We discuss the details of this hysteresis in 
Section 5. 

Fig. shows the pull-in and pull-out dynamics mechanism 
arising out of the gate dependent alteration of the two well 
potential. The pull-in occurs when a sufficiently high positive 
gate bias is applied between the gate and the source electrode 
until at a critical voltage Vpi, the local maximum between the 
two minima disappears. The pull-out requires a high negative 
gate voltage Vpo in order to break the Van der Waals adhesive 
bond. Described in terms of the potential landscape, the higher 
order flexural terms elevate the Van der Waals minimum until 
the latter becomes the metastable state and the cantilever 
configuration makes a reverse phase transition back to the 
horizontal coordinate. The resulting potential and the average 
angle are shown in that figure. Clearly the Vpo depends on 
the details of the Van der Waals potential and the higher order 
terms in the fiexural modes, and can be reduced by adjusting 
the contact geometry, material and surface properties such as 
its roughness and quality. 

IV. Sharp switching and breaking the Landauer 

LIMIT 

Let us try to work out the phase transition physics and the 
associated reduction in subthreshold swing. For a long stiff 
cantilever, the fluctuations about the mean can be ignored, so 
that P{0) oc becomes a delta function around the most 

probable angle set by the gate voltage Vg- We can expand 
around this extremum to give U « U{0*)-\-U"{A0)‘^/2, (each 
/ denoting an angular derivative with 0), so that the fiuctuations 
in the Boltzmann distribution look like a Gaussian, 

-U{e)lkBT ^ -U{e*)/kBT /2kBT (g) 

with standard deviation ^/kpTJU”. For a stiff cantilever 
compared to thermal fiuctuations, we can approximate this 
Gaussian as a delta function around AO = 0, which then 
integrated over angle gives us 

1=^ dEM{E)T(^E,e*{VG)) (3 Ft{E,Vg) (7) 


evaluated at a single gate-dependent angle 0*{Vg)- Differen¬ 
tiating with the gate voltage, we get 


dl 


2q 

h 



dEM{E) 


dT 


O Ft + T O 


dEp 


( 8 ) 


The tail of the thermal broadening function approaches the 
bandedge with an applied gate bias, making dEpIdVc ~ 
acP = aGq/kpT for an enhancement mode transis¬ 
tor, while the derivative of the first term in the bracket 
can be written as (dT/dVc) = {dT/d0*){d0'^/dVc)- 
The WKB tunneling term T = Tq exp [—gives 
dTjdO'' = —2KTdtair{0'^)/dO'^. From geometry, tair{0*) = 
h — LsinO* ^ L{0o — 6>*) for a long cantilever where 
0 < 6>* < 6>o ~ h/L 1. Then dtair{0^)/dO"" —F, and we 
get the normalized transconductance as the separable sum of 
electronic and conformational contributions 


9m 


d\nl 
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^conf 

i/m 


(9) 


assuming perfect gate control (ac = 1). The normalized 
transconductance is inversely proportional to the subthreshold 
swing to within a constant factor of InlO/g, so that the 
corresponding subthreshold swings add in parallel, reducing 
the overall value 


5-l = 5-l+Cn/ (10) 

It is easy to see from Eq. that the transconductance is 
already larger (and thus the subthreshold swing smaller) than 
the purely electrostatic limit /3, and that this extra contribution 
can diverge if either the cantilever is very long (k,L ^ 1), or 
if the cantilever can be made to collapse abruptly with gate 
voltage Vg so that dO"" /TVg '> 1. 

The physics of the evolution of 6>* is easily seen from the 
potential profile in Fig. The gate shifts the minimum of 
the first well (the metastable state) in an analog format, till 
the cantilever coordinate falls into the Van der Waals well. 
The minimum can be obtained by taking df/(6>*)/d6>* = 0. 
This gives us a transcendental equation graphically shown in 
Fig. as the intersection of two sets of curves, one of which 
shifts vertically with increasing Vg- For low values of Vg, we 
have three intersection points corresponding to three extrema 
- the two well bottoms and the bump in between. At high 


-Vg = O.OOV Vg = 2.50V - Vg = 5.00V 



Fig. 4. The movement of the minimum, is shown in the potential prohle. 
The Vq shifts the metastable local minimum of the first well, until the local 
minimum disappears and jumps into the global minimum of the Van der Walls 
well. 
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Fig. 5. The transcendental equation obtained by taking dU{0*)/dO* = 0 
is shown as the intersection of set of two curves. First set of curves (dashed 
red) are shifting vertically up as the Vq value is increased. For the low 
values of the Vq, there are three in-section points corresponding to the three 
extrema: the local minimum from cantilever, the global minimum from Van 
der Waals, and the local maximum between two wells. At a the high enough 
Vg, local minimum and the local maximum collapse at a point of inflexion 
with d?U{6*)/diO*^, followed thereafter by only a single, global minimum. 
This movement of the minimums is the main reason of the phase transition 
in the relays. 


enough Vg however, the wells merge and there is only one 
point of intersection to the right, the global minimum at the 
bottom of the Van der Waals well, and this is what drives 
the phase transition. The transition occurs when the potential 
valley corresponding to the first metastable state reaches a 
point of infiexion, corresponding to the dual condition 

dU{0^) 
dO^ 

_ 1 

d{0^Y 

We can simplify the algebra considerably by replacing 
the Van der Waals potential with a stiff parabola around its 
minimum, and ignoring the oxide thickness for the moment. 
We can then write 

u{e) « \ut^^e^ + ]^kvDw[o-eof 

egirAV^ 

L{eo-e) 2L{9o-9) ^ 

where kvDW = 57.14eI/^/cr^ around the Van Der Waals 
minimum. Setting its derivative equal to zero to extract the 
first metastable minimum, we get 

{UQ‘^^‘^ + kvDw)0*-kvDwOo = ^ 

which gives us the gate voltage dependence of the cantilever 
coordinate 0'^{Vg)- Setting the second derivative to zero as 
well to get the infiexion condition gives as an additional 
equation, and from the pair of equations we can then extract 
the destablization point 



= 0 

= 0 (11) 


rapidly grows away from the sticking configuration). However, 
near 0 ^ Oq to within a correction term AO ^ 21^11" ignored 
in our Gaussian approximation, we will get an additional pull 
in as k grows in strength. 

From the above equation, we can now extract the slope of 
the O'^iVc) curve, which enters the normalized transconduc¬ 
tance Qjyi 


gm = P + 


” 1 ” U^ef f^airA-L 

3f/e//(^*-05estab)(^O-^*) 


(15) 


where Ueff = + kyow- As explained earlier, we can 

make the transconductance really high by either ramping up 
the total dipole moment /i or by collapsing at the destabiliza¬ 
tion point of the cantilever, 0* ^ ^destab- Near that point, the 
equation simplifies 
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The second term can be interpreted as 

Qtot ^0 
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(16) 


(17) 


where Qtot = dUjdVG is the total charge stored in the 
cantilever. The normalized transconductance reaches infinity 
at that point, which is an artifact of the fact that we only 
included quadratic terms in our potential, and ignored thermal 
fluctuations around the expected minimum 6>* assuming an 
infinitely stiff cantilever. As we discuss in the next section, 
including thermal fluctuations brings in a finite temperature 
term oc pC averaged over all angles, which limits 

the subthreshold swing from its extreme zero value, although 
in practice this number can be (and has been measured to 
be) incredibly small. Eq. [T^ shows a second pull in near 
0"^ ^ Oq due to the Van Der Waals term, giving a second 
abrupt jump in transconductance. In Fig. we compare our 
results with and without the addition of VDW forces. In the 
absence of VDW, the performance is limited to 60 mV/dec 
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We get a destabilization near 0 ^ 0o/3 when the bending 
forces dominate (this happens when we’re far from the Van der 
Waals potential; the spring constant of the VDW term actually 
varies with angle if we take into account the nonlinearity that 


Fig. 6. (A Comparison of the potential profile with the addition of VDW 
forces (shown by blue curves) and without the VDW forces (shown by red 
curves). (B The angular movement of the cantilever with applied gate-bias 
with and without the VDW forces. (C) With the VDW forces the cantilever 
gets pulled down instantaneously shown in red, causing an improved SS of 
VlmV!dev. Without the addition of VDW forces from the gate, the fastest 
the voltage swing is limited to thermal limit of 60 mV/dev as shown in dashed 
black line. Thus proving the importance of the VDW forces at the nano-scale 
limit on the sub-threshold swing. 
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and the cantilever angle 0 increases slowly. With the inclusion 
of the VDW forces, the movement gets faster with the second 
potential minimum, lowering the SS to 12mV/dev. Once the 
gate voltage is large enough to remove the local maximum 
in Fig. the switching occurs with the tip of the cantilever 
getting pulled in and adhering with the drain electrode. 

As we scale the cantilever size, the gate to surface area 
reduces, making the capacitive forces weaker compared to the 
dipolar terms. For molecular cantilevers such as in sodium ion 
channels, the dipole moment dominates and is the main force 
behind the steep switching flM . We can see this transition 
from electrostatic to dipole driven correlated switching in 
Fig. [ 7 ] where we scale the dipole term with the cantilever 
length. 


V. Effects of finite temperature and einite scan 

RATE 


So far, we calculated the subthreshold swing at room 
temperature, which means that we needed to include thermal 
fluctuations especially around the inflexion point. Such fluc¬ 
tuations will average out the sharp transition over an angular 
width AO « y^2kBT/U" that we had ignored earlier in Eq. 
Near the inflexion point however, this approximation breaks 
down. Starting with Eq. invoking an equilibrium thermal 
distribution P for the configurational modes, we can redo the 
calculation to extract the normalized transconductance 


{,9m)conf — 



where the term in square brackets is the correlation term 
C{0,Vg)- For a long cantilever, dUldVc computes to 

dU _ llOPCairAVc 

dVc L{0o -0) ^ ^ 

Note that while 6>* depended sensitively on Vq, we are 
now dealing with an independent variable 0 that we will 
ultimately integrate over. The finite integration domain will 
keep the subthreshold swing finite, even near the infiexion 
point. However, the correlation term blows up at a particular 
angle, making the overall transconductance large but finite. 
The subthreshold swing has been measured to go down until 
as little as 0.1 mV/decade CD. 

After executing the voltage-gated phase transition, the can¬ 
tilever sticks to the drain. It is the pull in from this adhesion 
that makes the cantilever angle change abruptly and causes 
its subthreshold swing to plummet. However, this means that 
for the reverse cycle, we will need to tear the cantilever away 
from the drain by catapulting it with a high negative gate bias 
out of the deep VDW potential back to the weaker parabolic 
potential. This adhesion creates a strong hysteresis loop, which 
bodes ill for subthreshold swing, because the average slope 
over the loop becomes quite large even though the walls of 
the loop remain steep. 

Let us discuss how this hysteresis will arise out of our 
equations. The essential physics is governed by a voltage 
driven phase transition from a metastable to a stable state. 
For the forward swing, we move from a weak elastic potential 


Lever Movement 



Fig. 7. Comparison of the model with the addition and removal of the dipole 
moments.. With the additional of dipole moments to the model the switching 
characteristics change. Without the dipole moment the pull-in happens at 
the 2/3 of the gap, where as with the strong dipole moments, the transition 
becomes much more linear instead of parabolic and the sticking occurs at 1/3 
of the gap. The model captures the exactly the same switching behavior as in 
ESI when the dipole moment is not included and the same design parameters 
are used (shown as dashed.) 


to the Van Der Waals well, while for the reverse, we wait 
till the gate voltage pushes the Van Der Waals potential floor 
above the bottom of the elastic potential. It is clear that if the 
voltages varied very slowly near steady-state, i.e., given ample 
time, the cantilever would spontaneously stick by bending over 
and finding its global minimum, and thereafter the hysteresis 
would disappear. It is therefore essential to avoid a steady state 
calculation but invoke a finite scan rate Veit), and do our 
calculations assuming that the time constants of the cantilever 
are slower than the scan rate of the gate voltage. While this is 
trivially true in most experiments, in our model, this means we 
can no longer assume an equilibrium Boltzmann distribution 
P{0) (X but must solve the Fokker-Planck equation 

Eq.g with a time dependent gate voltage to extract the non¬ 
equilibrium behavior of the cantilever including its hysteresis. 
In Fig. we show the distribution of electrons under a slow 
scan-rate, where the cantilever configuration tracks the local 
minimum of the potential profile, and also the fast scan-rate, 
where the conformational distribution function lags behind 
the evolution of the potential, generating the aforementioned 
hysteresis. For a given scan rate, we let our numerical Fokker 
Planck evolve until a maximum time set by the scan rate, and 
sample the configuration and the current at that intermediate 
time to plot the nonequilibrium distribution function. 

The microscopic model of hysteresis can be generalized for 
any system that have three different time scales: fast, slow, and 
intermediate. Hysteresis occurs when the transition occurs in 
the intermediate time scale between the fast and the slow. In 
the case of the NEMFETs, the slow time scale is the relaxation 
speed of the cantilever coordinate and the fast time scale is the 
maximum possible value of scan-rate of the top-gate (giving 
electrons almost no time to relax or drift with the potential 
profile). When the top-gate is driven with a slower scan-rate 
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Fig. 8. The numerical solution of the Fokker-Plank equation coupled with 
the Landauer’s theory results in switching voltages that are dependent on 
the scan-rate. As the scan-rate (1/Tmax) increases, the switching voltages 
shifts to right, more positive values, causing a hysteresis. If the scan-rate is 
fast enough, the cantilever do not have time to follow the potential, resulting 
in a delayed version of the Boltzmann distribution. Tmax = 1 represents 
the case with no hysteresis where the scan-rate and the movement speed 
of the cantilever is comparable (Cantilever can follow the rate of change in 
the potential) (shown in A). Tmax = 0.25 has a faster scan-rate than the 
cantilever can follow resulting in a delay in transition causing it to switch at 
higher gate bias. Tmax is measured in units of scan-time and it normalized 
to the case with no hysteresis. 

than the speed of the cantilever response function, the adiabatic 
approximation works. However as the scan-rate speeds up to 
an intermediate value, the switching gets postponed more and 
more, causing a scan-rate dependent switching voltage shown 
in Fig. 

VI. Conclusions 

The hysteresis loop in the fabricated large devices (micron) 
and the smaller dimension biological and molecular devices 
(nanometer) vary and should not be ignored when scaling the 
devices to smaller sizes. Since the capacitive forces are weaker 
in the smaller devices, , the effect of the VDW dominates. 
As the devices get scaled aggressively to nanometer sizes, 
the nano-cantilevers will have wider hysteresis loops. This is 
where the dipole forces might help improve the performance 
of the devices, as shown in the Fig. [7] Also the smaller the 
cantilever length, the harder it is for it to bend, so that the 
hysteresis width gets bigger. Means to reduce the hysteresis 
loop include scaling the width geometry of the cantilever to 
favor the dipole forces or underplaying VdW forces, designing 
a pull back electronic circuit with multigated structures, ma¬ 
terials engineering to make the VdW well shallower, going to 
slower scan-rate or relying on ferroelectic phase transitions to 
liberate the cantilever from the Van Der Waals well |[30l . These 
effects however affect the subthreshold swing, so a proper 
optimization of the material and design landscape needs to 
be undertaken. The equations developed in this paper provide 
a starting point for such analysis. 

VII. 
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